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ABSTRACT 

Cooling rates are calculated for superfluid neutron stars of about one solar 
mass and 10 km radius, with magnetic fields from zero to about 10 1 4 gauss, 
when possible internal friction effects are neglected. Our results show that 
most old pulsars are so cold that thermal ionization of surface atoms would 
be negligible. At an age of a million years and with canonical magnetic fields 
of 10 12 gauss, the estimated stellar surface temperature is several thousand 
to a hundred thousand degrees. However, if we neglect magnetic fields and 
superfluid states of nucleons, the same surfaces would be about a million 
degrees. 
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COOLING OF PULSARS 

The hypothesis that a pulsar is a rotating neutron star now seems quite 
established. Almost all models which attempt to describe its observed slowing 
down or the mechanism which gives the observed radio pulses must ascribe to 
the neutron star a surface magnetic field in excess of 10 12 gauss. Except near 
their upper mass limit, the interior has been predicted to be mostly neutron 
superfluid as soon as the internal temperature becomes small next to a Mev 
(see, for instance, Ruderman 1968, Ginzburg 1969, and Baym'et.al. 1969). 

The effect of the magnetic field on the stellar surface region will be to reduce 
considerably the photon opacity. Superfluidity surpresses the normal heat 
capacity. Both effects, which have not been included in previous cooling rate 
calculations (Tsuruta and Cameron 1966a), greatly accelerate cooling rates. 

In this paper, emphasis is placed on consequences of these effects on the cooling 
of neutron stars during the stage when they can be observed as pulsars. 

For this purpose, we chose a neutron star model with the following 
properties; stellar mass M = 1.07 M 0 , radius R = 12.33 km, and central energy 1 
density pc = 7.39 x 10 14 gm/cm 3 . (These are medium mass models with the 
V y type nuclear potential of Levinger and Simmons and the composite equations 
of state with equilibrium compositions which were constructed by Tsuruta and 
Cameron (1966b) and subsequently used by Hartle and Thorne (1968) and other 
authors.) Recently substantial improvements were made on neutron star models, 
mainly by the use of nuclear potentials more realistic than those of Levinger and 
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Simmons. In these improved models, a neutron star of the same mass has 
slightly smaller size and higher densities than in our models. (For instance, 
in the most recent models by Baym, Pethick and Sutherland (1971), R = 7.8 ~ 10 km 
and p c =1~ 2.6 x 10 15 gm/cm 3 when M = 1.07 M 0 . (See also Ikeuchi, Nagata, 
Mizutani and Nakazawa 1971.)) However, the effect of these small differences 
on the cooling was found to be negligible compared with the effects of the 
magnetic fields and superfluidity, and the general conclusions of this report 
will not be changed by the use of more realistic neutron star models. 

The cooling of a neutron star generally proceeds through neutrino- 
antineutrino emission from the interior and electromagnetic radiation from 
the surface. In the earlier stages of the cooling history when temperatures 
are high neutrino emission dominates, while later at relatively lower tempera- 
tures photon emission is more significant. In our calculations we included 
neutrino emission from the plasmon process, the URCA process, the neutrino 
bremsstrahlung process, the photoneutrino process, the pair annihilation process, 
and the neutrino synchrotron process. We consider the history of a neutron star 
mainly after it has cooled down to several billion degrees, because profuse 
emission of neutrinos brings the star to that stage in a period much less than the 
age of the youngest pulsar. In the absence of magnetic fields and superfluidity, 
the URCA process is dominant, (the plasmon process competes with it at higher 
temperatures and bremsstrahlung at lower temperatures (Tsuruta and Cameron 
1966a)). Contributions by other mechanisms mentioned above turned out to be 
relatively minor. The major influence of the presence of magnetic fields in 
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this regime is to reduce the URCA neutrino luminosity in a neutron star. 

Nucleon superfluidity also suppresses the URCA process significantly (Itoh 1971). 
This is because the beta decay and electron capture of nucleons in the URCA 
process are hindered by the appearance of the superfluid energy gaps. (The 
other neutrino processes mentioned above are not affected by superfluidity, 
because only normal degenerate electrons, not superfluid nucleons, are in- 
volved.) Thus, for magnetic, superfluid neutron stars, the most dominant 
neutrino cooling mechanism is estimated to be the plasmon process and the 
total neutrino luminosity is somewhat reduced. 

Photon luminosity is described by the usual 

L = 4ttct R 2 T e 4 , (1) 

where a is Stefan's constant and T e is the effective (surface) temperature of 
the star. Convection is very unlikely in neutron stars (Tsuruta and Cameron 
1966a) , so that energy transfer will be governed by the opacity in the star. 

The total opacity K is generally expressed as: 

1 _ 1 1 

K “ K R + K c ’ (2) 

where K R and K c are the radiative and conductive opacity, respectively. 

Radiative opacity is due to various photoelectric processes (bound-free, free-free 
and bound-bound transitions) and scattering by electrons (Compton scattering 
at T > 5 x 10 8o K and Thomson scattering at lower temperatures). The 
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relative importance of these processes depends in a complicated way on density 
and temperature. Generally in the atmosphere where densities are lower, 
electron scattering gives the greatest contribution to the opacity at higher 
temperatures but photoelectric processes become most important at lower 
temperatures; in the inner degenerate regions electron conduction is the major 
contributor to the inverse opacity. 

In the presence of a superstrong magnetic field the opacity is greatly 
reduced (Canuto, Lodenquai, Ruderman and Tsuruta 1972, hereafter referred 
to as Paper II). To a very rough approximation the radiative opacity is altered 
by modifying the corresponding photon cross section according to: 

= a J°)’ for CO « co H , (3) 

where co is the radiation frequency, co H = eH/m e c is the electron cyclotron 
frequency, and ^(H) and <^(0) are the photon cross section with and without 
a magnetic field, respectively. This equation was derived in the case of Thomson 
scattering in earlier papers (Canuto 1970, and Canuto, Lodenquai and Ruderman 
1971). It is also appropriate for free-free transitions, which are essentially 
classical, and a very crude approximation for bound-free transitions when the 
magnetic field strength is about 10 12 ~ 10 13 gauss (see Paper II). Bound- 
bound transitions were found to be unimportant. We have applied the magnetic 
correction factor of Equation (3) to the zero field opacities previously obtained 
by Tsuruta and Cameron (1966a). The Rosseland mean was used to take account 
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of the frequency dependence. The radiative opacity, thus obtained, can be 
expressed as 

K r (H) = a R K r ( 0); with a R - (^j . (4) 

The conductive opacity is similarly expressed as 

K C (H) = a c K c (0) . < 5 ) 

The correction factor a c is generally a complicated function of a magnetic 
field, density and temperature (for instance, see Canuto and Chiu 1969). How- 
ever, for degenerate matter, the temperature dependence drops out and the 
factor a c decreases with increasing magnetic fields and decreasing densities. 

In the density regions where the conductive opacity plays a major role, the 
ratio K„(H)/K r (0) decreases typically by a factor of about several to a few 
hundred for a magnetic field strength of 10 12 to 10 13 gauss (see Paper II for 
details) . 

The total thermal energy of a neutron star with "normal" matter is derived 
from the small thermal tail of the Fermi distribution function of degenerate 
fermions, (neutrons, protons, electrons, muons and hyperons if they are present), 
and the kinetic energy of non-degenerate heavy ions. The energy of heavy ions, 
however, will not be important especially at temperatures significantly lower 
than ion lattice Debye temperature. For superfluid neutrons and protons, the 
specific heat is exponentially reduced. 
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c s = C n (T c ) 8.5 exp (-1.44 T c /T) ; for T « T c . (6) 

where C s and C n are the specific heat of the superfluid state and the normal 
state, respectively, and T c is the superfluid transition temperature. T c varies 
greatly for different assumptions of nuclear potential, effective mass density, 
etc. According to Takatsuka and Tamagaki (1971) the (neutron) energy gap may 
be between a few tenths and several Mev., so that the corresponding T c should 
lie between 10 9 °K and nearly 10 1 1 °K. If T c is as high as the upper limit, 
nucleon specific heat will be insignificant compared to that of the electrons. 

On the other hand, if we take the lower limit for the value of T c the contribution 
of nucleons will not be negligible until the internal temperature becomes less 
than about 3 x 10 8 °K. Both extremes are exhibited in our calculations to show 
the present measure of uncertainty in superfluid transition temperatures. 

Results are summarized in Figures. Internal temperatures and surface 
temperatures are plotted as a function of time in Figures 1 and 2, respectively, 
for varying strengths of magnetic fields. In the solid curves, neutron and proton 
thermal energies are neglected so that electrons are the major contributor to 
the total thermal energy. This corresponds to high temperature superfluidity 
(withT c ~ 10 11 °K). In the dashed curves, superfluidity is neglected and the 
heat capacity comes mainly from the thermal tail of the degenerate neutrons. 
Hence for a given magnetic field strength, the corresponding solid curve and 
dashed curve approximate the upper and lower limit, respectively, for the 
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heat capacity effects of superfluidity. The crosses indicate the points where 
photon luminosity takes over from neutrino luminosity as the main cooling 
mechanism. (Magnetic fields are expressed by the ratio 8 = H/H q , where 
H q = m e 2 c 3 /fie= 4.41 x 10 13 gauss.) 

In Figure 1, the internal temperatures are lower for solid curves (high 
temperature superfluidity) than for dashed curves (no superfluidity) at a given 
magnetic field strength (with H \o ). This is because the total thermal energy 
is decreased for superfluid matter. From Equation (4) we note that the ratio 
K r (H)/K r ( 0) decreases as (T/H) 2 and for sufficiently cool neutron stars the 
atmosphere becomes transparent for photons with certain directions and 
polarizations. When this point is reached the opacity is governed by electron 
conduction in the degenerate layers beneath the stellar surface. These layers also 
become increasingly transparent as the star becomes colder, because the 
degeneracy spreads to outer regions of lower densities with decreasing 
temperatures. 

It appears that during the earlier periods while the star cools by neutrino 
emission, surface temperatures are somewhat higher at a given age in the pres- 
ence of stronger magnetic fields. This is because the difference between the 
internal temperature and the surface temperature decreases due to the increased 
transparency of the strongly magnetized atmosphere, while the dominant neutrino 
luminosity is decreased. This effect becomes significant when the magnetic 
field strength becomes higher than about 10 13 gauss (see Figure 2). In the 



10 


presence of strong magnetic fields, the temperature drops fairly suddenly soon 
after electromagnetic radiation overtakes neutrino emission as the major cooling 
mechanism. This is because the atmosphere transparency increases at a 
growing rate as the temperature drops until the opacity almost vanishes. 

Exactly when this sudden cooling takes place depends critically on the field 
strength and the state of superfluidity. 

It is quite possible that the interior magnetic field is stronger than the 
value at the surface. The models indicated as (III) in the figures were con- 
structed to see this effect. In these models the magnetic field at the surface 
is H/H q = 0.025, the same as in models (II), but the magnetic field is increased 
by a factor of ten (times the surface value) beneath the surface where the con- 
duction predominates the inverse opacity. We see that the cooling rates are 
significantly increased for these models, more than the rates for the models 
(IV) witha constant magnetic field of H/H q = 0.1. 

An interesting result is that most pulsars may be very cold. Greenstein 
(1971) has pointed out that internal fricition may generate enough heat to main- 
tain the surface temperatures of certain pulsars somewhat above the values 
estimated here. (Such a possibility was first suggested by Cameron (1970) .) The 
magnitude of these frictional effects, however, depends upon the strength of the 
coupling between the charged and the superfluid components of the star, and 
quantitative estimates for all the possible ways this may occur are not yet 
available. The average age and field strength of pulsars appear to be a few 
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million years and about 10 12 gauss, respectively (see, for instance, Gunn and 
Ostriker 1970). For this age and field strength, the stellar surface tempera- 
ture is a little less than a thousand degrees for the solid curve (the exaggerated 
superfluidity) and about a factor of 10 3 higher for the dashed curve (without 
superfluidity), for models (II) in Figure 2. The actual value would be some- 
where between these two extremes. However, pulsars are certainly expected 
to have superfluid nucleons at internal temperatures much lower than the 
minimum calculated value of T c (~ 10 9 °K). It is estimated that the actual 
surface temperature is about a few thousand to nearly a hundred thousand de- 
grees (if frictional heating is negligible) . In the absence of a magnetic field and 
superfluidity, the surface temperature of a neutron star at the same age is 
almost a million degrees. 

For older pulsars, thermal ionization of atoms in the atmosphere may be 
negligible; an ionization energy of a few hundred eV is required at H ~ 10 12 
gauss (Cohen, Lodenquai and Ruderman 1970) while the atmospheric tempera- 
ture may drop to less than ~10 5 °K. (The heating effect of crust-core rotation 
slippage is the main uncertainty.) The neutron star surface is very likely iron, 
and this would make the ionization even more difficult. Therefore, if radio 
pulsars depend upon electrons streaming from the surface of a neutron star, 
decreasing surface temperatures below ^lO 5 °K may gradually turn off the 


pulsar (see Ruderman 1971) . 
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FIGURE CAPTIONS 

Figure 1. Internal temperature as a function of time, for varying strengths 

of magnetic fields. The models (I), (II), (IV) and (V) have the constant 
magnetic field of 6 = H/H q = 0, 0.025, 0.1 and 1, respectively. The 
models (HI) are for H/H q = 0.025 in the radiative opacity regions and 
H/H q = 0.25 in the conductive opacity regions. The solid curves 
approximate the upper limit for the heat capacity effects of super- 
fluidity. In the dashed curves superfluidity is neglected. 

Figure 2. Surface temperature as a function of time, for varying strengths of 
magnetic fields. The notations are the same as in Figure 1. 
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